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Abstract. Radiometeric CMB measurements need to be highly stable and this stability is best obtained with differential 
receivers. The residual 1// noise in the differential output is strongly dependent on the radiometer input offset which can 
be cancelled using various balancing strategies. In this paper we discuss a software method implemented in the Planck-LFI 
pseudo-correlation receivers which uses a tunable gain modulation factor, r, in the sky-load difference. Numerical simulations 
and experimental data show how proper tuning of the parameter r ensures a very stable differential output with knee frequencies 
of the order of few mHz. Various approaches to calculate r using the radiometer total power data are discussed with some 
examples relevant to Planck-LFI. Although the paper focuses on pseudo-correlation receivers and the examples are relative 
to Planck-LFI, the proposed method and its analysis is general and can be applied to a large class of differential radiometric 
receivers. 
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1. Introduction 

The dramatic progress achieved in the past decade in Cosmic 
Microwave Background (CMB) observations, particularly in 
anisotropy experiments, is strongly correlated with the remark- 
able improvements obtained in microwave and sub-millimetre 
detector technology, as well as in cryogenic t echnology (for a 
recent review see, e.g.. lBersanell i et ak 2002). Single-detector 
sensitivities of 0.1 - 0.4 mK Hz '^^ have been demonstrated 
for cryogenic operation of Indium Phosphide HEMT (High 
Electron Mobility Transistors) amplifiers (typically cooled at 
~ 20 K) in the range 25-100 GHz. At higher frequencies, sen- 
sitivities at a level NEP ~ 1 x 10 WxHz '^^ have been 
achi eved by spider web bolometers cooled to ~ 0.1 K (see, 
e.g., lLamarre..l997) . Moderate-size arrays of such detectors 
can today produce high resolution full-sky maps of the CMB 
with a high signal-to-noise ratio. These ultra-sensitive systems 
impose stability requirements that are proportionally stringent 
and call for highly optimised instrument design. In particular, 
the instrument needs to be immune at fiK level from the effect 
of parasitic signals introduced by non-idealities in the system, 
which would propagate as s ystematic errors in the final CMB 
maps jMennella et alJ2002h . 
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In coherent radiometeric systems one of the major concerns 
is the intrinsic instability due to gain and noise temperature 
fluctuations of the amplifiers themselves, typically well rep- 
resented by a 1 //-type noise spectrum. Differential receivers, 
such as the Dicke-switched scheme, reduce the impact of am- 
plifier instabilities in the measured signal with a fast (typically 
~ 100 Hz) switch between the sky input port and a stable ref- 
erence, sometimes given by another horn pointed at the sky. 
Dicke-type receivers have a long history in CMB observations 
and were successfully employed in the CO BE-DMR i nstru- 
ment that first detected CMB anisotropics JSmoot et alJll990l 
Il992h . 

In recent years, a scheme called "pseudo-correlation" ra- 
diometer has been introduced to improve over the classical 
Dicke scheme. As it will be discussed in detail, this design has 
a two-port front-end that allows a continuous comparison be- 
tween the sky signal Tsky and a stable reference signal Tret , im- 
proving the sensitivity by a factor V2 over a Dicke radiometer. 
In addition, fast (few kHz) phase switching provides immu- 
nity from back-end fluctuations. Different versions of pseudo- 
correlation designs are being used for the second and third 
generation of space-based radiometric instruments for CMB 
anisotropy: NASA's Wilkinson Microwave Anisotropy Probe 
(WMAP) and the Low Frequency Instrument (LFI) on board 
ESA's Planck mission. 
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In principle, a perfectly balanced pseudo-correlation ra- 
diometer is completely free from 1// effects. Residual sensi- 
tivity to 1// noise, as well as to other systematic effects, is pro- 
portional to the input offset AT^ff = T^^f - Tsky at the level of 
the first hybrid coupler in the radiometer front-end. In practice 
full balance is not achievable, so that a key instrument design 
objective is to minimise Aroff . 

In the case of the WMAP instrument, the radiometers di- 
rectly measure temperature differences betwee n sky signals 
from two widely separate d regions of the sky ilarosik et alJ 
l2003t iBennett et alJl2003l) . This is accomplished with a sym- 
metric back-to-back double-telescope system, and pairs of 
feeds which provide the two inputs to the pseudo-correlation 
front-end. In this scheme, the contribution to the offset from 
external signals is only of few mK (dominated by the CMB 
dipole of Galactic plane at low frequencies). The offset (of 
order AJoff <1 K) is dominated by second-order instrument 
asymmetries. 

The LFI radiometers, instead, measure differences between 
the sky, Tsky, and a stable internal cryogenic reference load 
(Tref ) cooled at about 4 K by the pre-cooling stage of the High 
Frequency Instrument (HFI) bolometer array in the Planck 
focal plane. This introduces an offset of order AFoff ~ 2 
to 3 K. The effects of this offset, however, are compensated 
by introducing a "gain modulation factor", r, which balances 
the output in the on-boa rd signal processing dBersanelli et al.l 
1 1995 ': Seiffert et al. 2002). Experimental re sults from advanced 

LFI prototypes (Meinhold et al 1998: Tuovinen et aL20 00). as 

well as analytical calculations, show that great immunity from 
1 // effects can be obtained when the value of r is accurately 
selected. 

Two different approaches can be used to select the value 
of r. it can either be controlled in hardware by adjusting (in 
principle in real-time) a variable gain to achieve the null-output 
condition; or it can be set as a controllable software parameter, 
giving more flexibility to the system. The latter approach was 
adopted for the LFI. 

As we shall see in more detail, the best estimate of r can 
be obtained based on the radiometer data themselves, collected 
in a raw, undifferenced form. In a space application this may 
impose non-trivial demands on the telemetry rate. If the data 
available to calculate r are limited by telemetry constraints, this 
may limit the accuracy of the determination of r and therefore 
limit the instrument stability. 

A detailed analytical study of the impact of non-idealities 
in the LFI radiometer, in partic ular 1 / f noise effec ts, has been 
carried out in a previous work jSeiffert et alJ2002l) . In this pa- 
per we discuss in detail the required and obtainable accuracy of 
r, and its effect on the measurement quality. Although we fo- 
cus on offset balancing in pseudo-correlation radiometers and 
discuss some examples in the context of Planck-LFI, the con- 
cepts and formalism discussed in our paper are general and 
appUcable to any switched radiometer like unbalanced Dicke 
receivers. 

After a short description of the radiometer concept, we dis- 
cuss the issue of offset balancing and derive the optimal r for 
1// noise minimisation (Sect. |3}, demonstrating its effective- 
ness by applying it to representative simulated data streams. In 



addition we discuss the impact of the choice of r on other sys- 
tematic effects and derive the required accuracy on the calcula- 
tion of r. In Sect.l^we present and compare various approaches 
for calculating r and show how these can be used in the con- 
text of space experiments; an example of the application of this 
offset-balancing technique to laboratory radiometer data is pre- 
sented and discussed in Sect|5] Finally, in Sect.|6lwe discuss the 
impact of systematic effects of instrumental and astrophysical 
origin on the calculation of r. 

2. Pseudo-correlation differential radiometers: 
basic concepts 

The Dicke switched radiometer represents a typical implemen- 
tation of a differential receiver. In its simplest form a switch 
located in the front-end commutes rapidly between the sky and 
the reference horns so that a sequence of sky-load signals is de- 
tected and differenced at the output of the amplification chain; 
the radiometer sensitivity per unit integration time is given by: 

Ar™, = 2%, (1) 

where Tsy, = T^i^y + T^oi^g is the system temperature, r„oise 
represents the radiometer noise temperature and /3 is the band- 
width. 

If the switching frequency is sufficiently high so that the 
gain can be considered nearly constant in a sky-load cycle, then 
the susceptibility to gain fluctuations is very much improved: 

Arcain _ AG Tsky - T^ref 
sys ^ i sys 

where Arcain represents the rms signal variation caused by gain 
fluctuations and Tsky, T^et represent the sky and reference load 
temperatures, respectively. 

Although this scheme is effective in reducing 1// noise in 
the final measurements, the presence of a lossy active compo- 
nent in the front-end increases the radiometer noise and may 
introduce additional 1// components. These limitations, how- 
ever, can be overcome with a modified version of the Dicke 
scheme, the pseudo-correlation radiometer, that has been re- 
cently adopted for the WMAP and Planck-LFI instruments. 

In Fig. n we show a schematic of a pseudo-correlation ra- 
diometer in its simplest form. The sky and reference signals are 
summed by a front-end 180° hybrid, amplified by two parallel 
amplification chains and then separated by a second hybrid, de- 
tected and differenced. In this way the 1// noise from the RF 
amplification stages is the same in both signals, so that they can 
be eliminated at first order by differencing. 

To reduce the effect of instabilities in the back-end elec- 
tronics a fast switching between the two inputs is implemented 
in the two radiometer legs; as an example of application of the 
switched pseudo-correlation scheme, Fig.|2shows a schematic 
of the radiometers used in the Planck-LFI instrument. 

In the front-end part (see top part of figure) the radiation en- 
tering the feed-horn is separated by an OrthoMode Transducer 
(OMT) into two perpendicular linearly polarised components 
that propagate independently through two parallel radiometers. 
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Fig. 1. A pseudo-correlation radiometer in its simplest form. 
This scheme allows a differential measurement without requir- 
ing an active switch before the first RF gain stage. 



In each radiometer, the sky signal and the signal from a sta- 
ble reference load at ~4 K are coupled to cryogenic low-noise 
HEMT amplifiers via a 180° hybrid. One of the two signals 
then runs through a switch that applies a phase shift which os- 
cillates between and 180° at a frequency of 4096 Hz' (the 
second phase switch is present for symmetry on the second ra- 
diometer leg but it does not introduce any phase shift). The 
signals are then recombined by a second 180° hybrid coupler, 
producing a sequence of sky-load signals at the output alternat- 
ing at twice the frequency of the phase switch. 

In the back-end of each radiometer (see bottom part of 
Fig. 13 the RF signals are further amplified, filtered by a low- 
pass filter and then detected. After detection the sky and refer- 
ence load signals are integrated and digitised before sending to 
ground. Note that the current Planck-LFI available telemetry 
allows the downloading of the total power data streams for all 
the detectors, so that the sky-load differecing will be performed 
"off-line" during ground data analysis. 



3. Offset balancing and 1// noise minimisation 

In differential radiometers it is desirable to maintain the offset 
between the sky and reference signals as small as possible in 
order to minimise the effect of gain instabilities (see Eq. (|2j)- 
An offset introduced by an internal reference load can be bal- 
anced before differencing either by a variable back-end gain 
stage with a feed-back scheme to maintain the power output 
as close as possible to zero, or by by multiplying in software 
one of the two signals by a so-called gain modulation factor. 
In both cases the differential radiometer output can be written 
as: 



Pout = aGtotksP [fsky + Tn- r (f ref + Tn)] , 
^sky — ^sky/^sky + (1 ~ ^sky^"^phys' 
fref - Tref/Lref + (1 - i'ref)^phys' 



(3) 



where Lsky and Lref are the insertion losses of the front-end sky 
and reference load antennas, Tphys is the physical temperature 
of the receiver front-end, a is the detector proportionality con- 
stant, Gtot the radiometer total gain, the Boltzmann constant. 



' This frequency corresponds to a full 0-180° cycle, so that the rate 
between two subsequent switch states is 8192 Hz 



20 K Front end 

Ref. load 
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Fig. 2. Baseline LFI pseudo-correlation radiometer One of the 
two 180° phase switches in the front end switches rapidly the 
phase of the propagating signal, thus producing a sequence of 
sky-load signals at the output of the second hybrid. In the warm 
back-end the signal are further amplified, detected and digitised 
before telemetry. 

/3 the radiometer bandwidth and r the gain modulation factor. 



From Eq. Q it follows that po 

Tsky + Tn 



for r - where 



(4) 



In the next section we will show how the above condition 
of zero power output represents a very good approximation of 
a balanced radiometer and is very effective in reducing the ra- 
diometer susceptibility to gain fluctuations to very low levels. 

In the case of Planck-LFI the radiometer offset has been 
balanced through a software scheme that uses the undifferenced 
total power data to calculate r, which is expected to be constant 
on timescales as long as several days. This scheme has been 
preferred to a hardware gain modulation because it simplifies 
the radiometer hardware and avoids potential systematic errors 
from the variable gain stage. 

3.1. Analytical derivation of the balancing condition 

In this section we calculate the value of the gain modulation 
factor that cancels the radiometer susceptibility to front-end 
gain fluctuations and show that this value is very well approxi- 
mated by the the value in Eq. (0}. 

Before proceeding into the details of the calculation we 
briefly examine the expected magnitude of gain and noise tem- 
perature fluctuations in the total power noise streams; more 
details concerning the analytical trea tment of LFI pseud o- 
correlation radiometers can be found in ISeifl^ertet alJll2002l) . 

3.1 .1 . Noise fluctuations in total power streams 

Cryogenic HEMT amplifiers are knowi i to have 1/f fluc- 
tuations in gain and noise temperature (|Pospieszalskilll989l: 
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IWollacklfl99l ljarosiklll99dl) . The level of these fluctuations 
can vary considerably among amplifiers and depends on the de- 
tails of device fabrication, device size, circuit design, and other 
factors. Because of this, we adopt an empirical model for the 
fluctuations also described in Seiffert et al. (2002). 

In particular we write the 1// spectrum of the gain fluctua- 
tions as 



AG(/) ^ C_ 
G 



(5) 



where 0.5 < a < \ and C represents a constant normalization 
factor. Similarly, we can write the noise temperature fluctua- 
tions as 



Arn(/) 



A 



(6) 



where A ~ C/(2A^s) is the normalization constant for noise tem- 
perature fluctuations. In the case of the Planck-LFI radiometers 
we use the estimates A ^ 1.8x10"^ for the 30 and 44 GHz and 
A ^ 2.5 X 10"^ for the 70 and 100 GHz radiometers. 

3.1 .2. Derivation of tine optimal gain modulation factor 
for 1// minimisation 

Following the approach described in lSeiff'ert et al ] (l2002h it is 
possible to derive analytical formulas for the knee frequency 
relative to 1// fluctuations of the differenced data streams in- 
duced by front-end gain and noise temperature fluctuations. 

If we assume that gain and noise temperature fluctuations 
are correlated then we obtain the following relationship for the 
knee frequency: 



- r(fref + (1+A/C)r„)]} 
which is zero for the following value of r: 

fsky + (l+A/C)rn 
fref + (1 +A/C)T,' 



(7) 



(8) 



If we relax the hypothesis of completely correlated gain and 
noise temperature fluctuations then we find: 



C 



(1 - rfTl(A/Cf] 



((f sky + T^n) - Kfref + ^n)) 
1/2] I/O- 



(9) 



which is minimised (but not nulled) by the following value of 
r: 

^ , _ (f sky + T„)(T,,f + Tn) + (T,A/Cf 

' (f,ef + r„)2 + (r„A/c)2 ■ ^ 

Note that in both Eq. (|8j and ilOi the optimal value of r is 
dependent on the ratio C/A ^ 2 y/Ns where A^s is the number of 
amplifier stages; if C/A » 1 then we have that in both cases 
the optimal value is approximated by (given by Eq. @). 



In Fig.|3we show a plot of the knee frequency calculated 
using Eq. and (|9} versus the relative accuracy - (f - 
r*)/r* for parameters typical of the 70 GHz radiometer chain 
(which drives the accuracy requirements in the calculation of 
r). 

Two main consideration may be drawn from the results 
in Fig. 13 first, if r = then the radiometer knee frequency 
is not nulled also if we assume complete correlation between 
gain and noise temp erature fluctuations. In fact it can be shown 
dSeiffert et alJ2002i) that gain fluctuations are cancelled at first 
order and the residual 1// noise comes from noise tempera- 
ture instabilities. Second, if r,*, is known with an accuracy of 
the order of ~ 1.5% then the resulting knee frequency is be- 
low the design value (50 mHz in this case) independently from 
the hypothesis of correlated or uncorrelated fluctuations in the 
front-end amplifiers. 




-0.03 -0.02 -D.Ol 



O.Dl 0.02 0.03 



Fig. 3. Radiometer knee frequency versus the relative accuracy 
6r/r - if - r^^jr^ for the 70 GHz Planck-LFI channel. The 
solid and dashed curves represent the values of the knee fre- 
quency calculated assuming correlated and uncorrelated fluc- 
tuations in the front-end amplifiers. The dashed-dotted 50 mHz 
line indicates the Planck-LFI design value for the knee fre- 
quency, while the shaded area represents the range of accuracy 
that is compliant with this design value. 



In Fig. 0] we show a simulation of the total power and of 
the differenced noise stream for a 30 GHz LFI radiometer in 
the hypothesis of correlated gain and noise temperature fluctu- 
ations. The simulation has been performed considering typical 
values for system noise temperature and front-end losses, ~ 1 K 
contribution from the telescope (emissivity of 1 % for each re- 
flector at 50 K physical temperature), a physical temperature of 
4.8 K for the reference load and a physical temperature of 20 K 
for the LFI focal plane. 

In the upper graphs of Fig0]we show the sky and reference 
total power data streams, while in the middle graphs we com- 
pare the differenced noise streams obtained with r - \ (left 
panel) and r - r* (right panel); the graphs show that in both 
cases most of the 1 // instability is removed, although in the 
r - \ case the sky-reference offset is retained. A comparison 
between the two amplitude spectra (lower panels) reveals that 
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Fig. 4. Simulated noise streams with parameters representative 
of a a 30 GHz LFI radiometer The top graphs show the total 
power data (sky, left panel, and reference load, right panel), the 
middle graphs the differenced noise with r = 1 (left panel) and 
r - and the bottom graphs show the corresponding ampli- 
tude spectra. 



in the radiometer couple with the measured output essentially 
through the following mechanisms: 

- variations in thermal noise induced by resistive front-end 
passive components (feed-horn, OMT and first hybrid); 

- amplifier gain oscillations (in the front-end and in the back- 
end); 

- amplifier noise temperature oscillations (in the front-end^). 

The radiometer susceptibility can be written in terms of a trans- 
fer function, (p, that links a physical temperature variation, 
(ST^phys, to the corresponding variation in the measured output, 
^rout, i.e.: 6Toai - (f> x dTpi^ys. This func tion has beeri calcu- 
lated following the approach discussed in lSeiffert et al 1 (l2002h 
assuming complete correlation between the temperature fluctu- 
ations at the level of the various radiometer components. This 
assumption is a reasonable one in the context of the Planck- 
LFI radiometers where physical temperature variations induced 
in the 20 K and 300 K stages will be dominated by relatively 
"slow" components with a very small contribution from har- 
monics with a frequency higher than ~5 mHz. 

Under these assumptions the susceptibility to temperature 
variations in the 20 K stage (</ife) and in the 300 K back-end 
stage (<^be) can be written as: 

#E = isky {(1 - i^ky) - Kl - L;^!) + 

In(lO) r- - 1 SG^E 

+ —nr- Tsky + T^^ - r{T,^f + T^,^) — + 

+ (l-O^j, (11) 
phys J 



if the ideal gain modulation factor is applied then the final knee 
frequency is much lower compared to the r = 1 case^. 

In Tabled we provide representative estimates of for all 
the frequency channels of Planck-LFI radiometers. 



30 GHz 


44 GHz 


70 GHz 


0.936 


0.953 


0.971 



Table 1. Estimates of r,t for Planck-LFI radiometers. 



3.2. Susceptibility to other systematic effects 

Apart from 1/f noise, other instrumental systematic effects can 
be mitigated by balancing the offset in the radiometer output. 
In Planck-LFI, for example, systematic effects can be expected 
from thermal variations of the 20 K and 300 K stages and from 
input bias fluctuations of the front-end amplifiers. 

Let us consider, for example, the effect caused by thermal 
instabilities in the front-end and in the back-end stages of the 
Planck-LFI radiometers. Variations in the physical temperature 



Note that in the case r = the knee frequency could not be re- 
solved due to data stream length (1000 s) that limits the frequency 
resolution to ~ 5 mHz 



= L 



sky" 



In(lO) 



da' 



,dB 



10 yOTphys dTphys) 
X [f sky + Tsys - Kfi-ef + Tsys)] , (12) 

where Gpg, and are the gain in dB and the noise temper- 
ature of the front-end amplifiers, Ggg the gain of the back-end 
amplifiers, a*^^ the diode constant and Tphys the physical tem- 
perature. 

Note from Eq. il2\ that the susceptibility to back end tem- 
perature fluctuations is cancelled if the radiometer output is 
balanced (i.e. ^y^^ = for r = r^); this is not only true for 
temperature variations, but for also for any instrumental effect 
causing instability in the back-end amplifiers. 

The susceptibility to front-end temperature variations, in- 
stead, is cancelled for a value of r which is slightly different 
from (its expression is not reported here, but it can be de- 
rived from 0FE = in Eq. (^3)- In general the effect of any 
given systematic effect in the radiometer can be cancelled at 
first order using a particular value of the gain modulation fac- 
tor; in practice, however, the choice of the "best" value of r to 
use in the difference is driven by the most critical expected ef- 
fect and by the feasibility to calculate r from the measured data 
with the required accuracy. 



The ~35 dB gain in the front-end makes negligible the effect from 
fluctuations in the noise temperature of the back-end amplifiers 
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In Planck-LFI, for example, with the condition r = it is 
possible to minimise 1// amplifier noise and effects from the 
back-end while keeping the radiometric susceptibility to other 
systematic effects at a level which allows to define realistic re- 
quirements on the thermal and electrical stability at the inter- 
faces between the satellite and the instrument. 

3.3. Required accuracy in the calculation ofr*^ 

The main goal for balancing the radiometer output is to min- 
imise the residual 1// noise; therefore it is natural that the re- 
quirement on the gain modulation factor accuracy must be de- 
rived from the maximum knee frequency allowed in the exper- 
iment. 

From the curves shown in Fig.l^we have derived a require- 
ment on the accuracy of the gain modulation factor of the order 
of - ±1.5%. 

In some cases, however, the accuracy needed in the cal- 
culation of the gain modulation factor may be determined by 
the level of other effects; in Planck-LFI radiometers, for ex- 
ample, the value r,*, will not only minimise 1// noise, but also 
effects from back-end thermal instabilities; in this case the need 
for a very low radiometer susceptibility to back-end tempera- 
ture fluctuations imposes a tight requirement on the accuracy 
of the gain modulation factor, ranging from ~ ±1% at 30 GHz 
to ~ ±0.3% at 70 GHz. 

4. Calculation of r* from measured data 

During a long-duration CMB measurement the radiometer 
noise properties are subject to slight, slow variations in time 
due to thermal/electrical drifts and ageing of the radiometer 
electronic components. If the offset balancing is done in soft- 
ware, then the value of can be recalculated and updated (if 
necessary) in the data reduction software to compensate for 
such changes. 

As we show in detail in this section, the gain modulation 
factor can be calculated from radiometer data acquired in "total 
power mode", i.e. before differencing is performed. Although 
the current telemetry bandwidth available to the LFI instrument 
allows to download the complete total power data streams, it 
is possible to calculate r from a limited portion (about ~15 
min) of these data. This implies that it could be possible to 
perform the sky-load differencing on-board (with a saving of 
~50% in telemetry bandwidth) providing margin for possible 
non-nominal operation scenarios. 

In this section we discuss the following three calculation 
methods and compare them from the point of view of the ob- 
tainable accuracy: 

1 . calculated from the ratio of the average sky and reference 
load levels; 

2. r,* calculated from the ratio of the sky and reference load 
standard deviations; 

3. calculated by minimising of the final differenced data 
stream knee frequency. 

In the following discussion we will consider only the presence 
of white+1// noise in the radiometer data streams and assume 



the availability of 15 min of total power data (as in the case of 
Planck-LFI) in order to evaluate the efficiency of each method. 
The effect of the presence of the astrophysical signal and of 
other instrumental systematic effects is discussed in Sect.|6| 

4.1. Calculation from average signal level 

The simplest way to calculate r from the total power data is to 
use directly the definition provided by Eq. 0, i.e. to take the 
ratio between the average signals acquired when looking at the 
sky, Vsky,; - G,(7'sky + T^j), and when looking at the reference 
load, Vi-ef,/ = Giifiei + Tn,i)- In this case r is given by: 



Hill ^lef.i 



(13) 



where represents the number of samples available in the data 
stream. Now we know that although gain and the noise tem- 
perature will display l/f fluctuations, thanks to the pseudo- 
correlation scheme they will be the same (at first order) in both 
the sky and reference signals. Therefore for each sample, /, we 
can write: G, = G + 6G^.^^ and T„i - T„ + ST^^f, where G 
and T„ represent the average gain and noise temperature level. 
Furthermore there will be a white nose component ST'^^^^^^ 
such that {6T^^^^.^^} = and cr^^^^^ = (f^j^^^^.^^) + r„)/ y/pT. 

Let us now evaluate the impact of gain and noise temper- 
ature fluctuations separately. If the total power noise streams 
contain l/f gain fluctuations then we can write the sky and 
load data streams as: 



Vsky(ref),/ = (G + 5G,'^^)(f sky(,-ef) + T„) + 6T^ 



sky(ref) ' 



therefore the sums in Eq. J13> have the form: 

N 

^ V'sky(ref),i = GA^ skydef) + T'n) X 



(14) 



(15) 



N ^ 



nWN 



N 

Zj ^ky(ref)- 
;=1 



Similarly, in the case of l/f noise temperature fluctuations 
then Eq. (I15> becomes: 



^ Vsky(ief),/ - GA(fsky(i-ef) + r,,) : 



1 + A" 



ST"f 



sky{i'ef) 



(16) 



N 

ST. 



1=1 



sky(ref) " 



In both Eqs. il5\ and J16t the two sums on the right-hand 
side tend to zero for large values of A, so that the ratio of the av- 
erage signals — » rj^ for A — > oo, which implies that Eq. ( I13t can 
be used also in presence of l/f noise provided that the statistics 
is sufficiently large. 

In Fig.|5lwe show the result of a numerical simulation of the 
convergence of Eq. il3\ to the theoretical value r^. Each curve 
represents the value of Sr/r^ versus time, averaged over 50 
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equivalent noise realisations of the 30 GHz Planck-LFI chan- 
nel. The different curves are relative to different levels of 1// 
noise in both gain and noise temperature (A = 0, i.e. white 
noise only, A - 10"^, A - 10 A - 10"^ i.e. low, interme- 
diate and high level of 1// noise). The figure shows that the 
presence of 1// fluctuations determines a slower convergence 
with respect to the white noise case. Even in the worst case, 
however, (which represents an unrealistic case with respect to 
the expected LFI performances) the accuracy obtained with 15 
minutes of data is much better than 0.01%. 



1 .ooE-oa 



1 .OOE-04 



1 .OOE-05^ 



1 .00E-0& 



1 .OOE-07 




10 100 
Time (s) 



1000 



Fig. 5. Effect of gain and noise temperature 1// fluctuations on 
the calculation of r from the ratio of the average signal levels. 
Each curve is an average of 50 noise realisations relative to the 
30 GHz Planck-LFI channel. 

In order to calculate properly the ratio r from Eq il3\ . it is 
necessary that any "zero-level error" introduced by the back- 
end electronics is small. Let us evaluate the effect of a constant 
offset SV added to both states in the switching output of the 
radiometer; the voltage measured in the two switch states can 
be written as: 



sky(ref) 



^ G(f sky(ref) + r„) + 6V. 



(17) 



If the offset is unknown, it introduces an uncertainty given 



by: 



6r : 



Vsky 



Vsky + 6V Vsky - Vi-ef 



Vref + SV Kef(Kef + SV) 



-6V. 



(18) 



From Eq. J18> it is straightforward to obtain the maxi- 
mum values of the unknown offsets that would produce an 
uncertainty on r equal to the required accuracy. In the case 
of Planck-LFI these values (in temperature units) range from 
~ ±0.5 K at 30 GHz to ~ +2.34 K at 100 GHz, which are 
compatible with the current electronics design. 

Finally we underline that the above treatment is valid in 
the hypothesis of detector linearity and negligible gain com- 
pression. In general, however, these conditions are largely met 



in CMB anisotropy experiments, where the signals to be mea- 
sured are small. In the case of the Planck-LFI radiometers, 
for example, the radiometer output power is of the order of 
-25 dBm, which is much smaller than the -1-dB compres- 
sion point of the HEMT amplifiers that is of the order of 
~ -15 dBm and falls well in between the range of detector 
Hnearity (~ -35 H- -15 dBm). 

4.2. Calculation from standard deviations 

A second approach is to calculate r from the ratio of the stan- 
dard deviations of the total power noise streams, i.e.: 



ref 



Sill (Vsky,,- - <Vsky>)' 



Zii (Kef,,- - (Vrf))' 



1/2 



(19) 



An advantage of this method compared to the previous one 
is that the ratio of standard deviations is insensitive to spurious 
offsets introduced by the back-end electronics; let us evaluate 
now the accuracy that is obtainable in presence of noise. In the 
case the data streams are characterised by white noise only then 
we have that cr^y(,.gf^ ^ (fskycref) + T„)/ ^J^t for ^ oo, so that 
for large values of N Eq. M9\ approximates r*^. 

Of couse, it is also well known that the standard deviation 
of an infinite time stream of l//-type noise formally diverges. 
For the case under consideration here, the ratio of standard de- 
viations of finite sets of white noise plus some l//-type noise, 
we examine below the details of accuracy and convergence. 

Let us consider first the effect of gain fluctuations; from 
Eq. (I15> we have that 



sky(ref),/ ~ 



/=i 

- G(fsky(ref) + ?"n)^l 
+ <^^s*y'(i-ef))A'' 



SG'lf 



(20) 



and, consequently. 



(Vsky(i-ef),/ ~ (Vsky(i-ef))) ~ 

1=1 

N 

(f skyCreO + T,f ^ [SG'/ ^ - {6G"^)f + 
1=1 

N 

l,"^ sky(rel),/ skyCietyy * 



(21) 



(A-1) 



Considering that cr^J^^^.^f^ ^ (f sky(ref) + To)/ ^J^t then it 



follows that 



" sky(i-ef) 



(y'skyCref) + To) 



// + 



(22) 



so that the ratio of the total power standard deviations con- 
verges to for large A^. 
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Let us now consider the presence of 1// fluctuations in the 
noise temperature. In this case Eq. illi has the form: 



(^sky(ref),i " (^sky(ref))) = 

N 

G'f,{ST:j/-(ST!Jf))\ 

;=1 

2lj \" sky(ref),! N"^ skyCref)// 
/ WN Y 



2 N- 



(23) 



i=l 

(N-l) 



It follows that the ratio of standard deviations converges, in 
this case, to: 



1// ^ '-'^ sky 



(24) 



This clearly shows that in presence of noise temperature 
1// fluctuations the ratio of standard deviations does not con- 
verge to Tq. In particular with increasing amplitude of the 1// 
noise temperature fluctuations (i.e. for increasing values of A) 
we have that Gcri/^ » cr*'^, „ and the ratio of standard devia- 

i n sky(ref) 

tions tends to 1 . We can therefore calculate the following lower 
limit for the relative accuracy on r obtained by this method: 



Sr A- 



1 



(25) 



In Fig.|6lwe show the results of numerical simulations rel- 
ative to the 30 GHz Planck-LFI channel showing the conver- 
gence of the ratio of noise standard deviations (Ea.ll9> for var- 
ious levels of 1// noise. The two dashed curves represent the 
limits for A — > (i.e. white noise only) and A — > oo. In the first 
case (white noise), because the standard deviation of random 
samples tends to the standard deviation of the gaussian distri- 
bution as dcr/a = 1 / "^V, we have that the accuracy calculated 
by the ratio of the standard deviations increases with the num- 
ber of samples as 6r/r - 2/ y/N. In the second case we have 
that the ratio of standard deviations does not converge to but 

l-r* 

to The figure clearly shows that for increasing values of A 

''o 

the behaviour changes between these two limits. 

When the 1// noise level is high it is still possible to esti- 
mate r from the ratio of the white noise rms that can be derived 
either by fitting the high frequency end of the noise power spec- 
trum or using a proxy of the white noise level by means of the 
two-point variance furiction cr^{i-j) = {(xi-xj)^} as described 
bv l.Tanssen et alJ lll996l) . In this case r will converge to rj*, with 
an upper limit accuracy of ~ 2/ Va^ (see Fig.|6j. 

4.3. Calculation by minimising the final knee frequency 

A third approach that we consider here is the optimisation of 
the gain modulation factor by defining a window [rmin, rmax] 
around r* and then finding the "best" value of r that minimises 
knee frequency of the differenced noise stream. 



0.1 



0.001 




10 100 
Time (s) 



Fig. 6. Efifect of noise temperature 1// fluctuations on the cal- 
culation of r from the ratio of standard deviations signal levels. 
Each curve is an average of 50 noise realisations relative to the 
30 GHz Planck-LFI channel. 

Although this method probes directly the 1/f noise charac- 
teristics of the final differenced data, there are some limitations 
and caveats that are worth mentioning: 

1. frequency resolution and accuracy. If Af represents the 
length (in time) of the data available for analysis, then 
the absolute minimum frequency that can be resolved in 
Fourier space is /min ~ 2/Af. For Af ~ 900 s (15 minutes) 
we have that /min ~ 2 mHz. Although the data length can 
be enough to resolve frequencies ~ 1 mHz, in order to de- 
termine /k we must be able to recognise the frequency at 
which the 1/f noise has an amplitude equal to the rms value 
of the white noise. Because noise dominates, each individ- 
ual sample has 100% uncertainty so that we must rebin the 
data with a consequent loss in frequency resolution. 

2. Presence of astrophysical signal in differenced noise 
stream. When the sky and load data streams are differenced 
taking into account a value of r in the interval [rmin, '"max], 
the differenced noise stream will also contain the astrophys- 
ical signal (dipole, galaxy, etc.) at all the harmonics of the 
scan frequency (~16 mHz for Planck). This implies that 
with the real data the calculation of the knee frequency 
will be affected, in principle, by the power at the first 10 
harmonics or so of the scan frequency. Therefore some as- 
sumptions on the expected signal are probably necessary to 
remove it from the differenced data before computing /k. 

5. Application to laboratory radiometer data 

The gain modulation strategy described in the previous sections 
is routinely applied with success in the Planck-LFI prototype 
radiometers. Here we present some results obtained with exper- 
imental data from the 30 GHz LFI Elegant BreadBoard (EBB) 
radiometer assembled at the Jodrell Bank Observatory (JBO) 
in the framework of the Planck collaboration. Further details 
about the EBB testing campaign of the Planck-LFI radiometers 
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will be discussed in a dedicated set of papers which is currently 
in preparation. 

The EBB 30 GHz radiometer is comprised of a front- 
end module (FEM) manufactured at JBO and a back-end 
module (BEM) manufactured at the University of Cantabria 
(Santander, Spain) connected by waveguides approximately 
0.75 m long. Each FEM channel includes a HEMT phase 
switch which has also been designed and manufactured by 
JBO. In the tests described here the FEM was operated at 20 K 
while the BEM was maintained at room temperature. 

Because the aim of these tests was to verify the validity of 
the differential concept rather than to achieve flight-like perfor- 
mances, the test setup was kept as simple as possible; in partic- 
ular the front-end phase switches were operated at 280 Hz (in- 
stead of 4096 Hz as in the LFI baseline) and the 300 K thermal 
environment was not maintained at a high degree of stability. 

Before the stability tests were carried out, the BEM white 
noise floor was measured to be significantly less than the 
noise level with the FEM and BEM connected together Fifteen 
minute tests were then carried out with one of the phase 
switches kept in a fixed state. A 280 Hz square wave switch 
waveform for the other phase switch was generated by the data 
acquisition card counter outputs, under Lab VIEW control. The 
detected signals were therefore 280 Hz square waves, which 
were sampled at 504000 samples per second. After discarding 
20% of the data (10% from each edge of the square wave), each 
switched state was averaged to obtain a single number, Ai, Bi, 
A2 or B2, where A represents the 'on' state for each channel 
and B the 'off' state. 

In the left graph of Fig. we show the total power levels of 
one of the two output radiometer detectors, the top trace repre- 
senting the A states and the lower the B states. The small inset 
shows how the actual averaged diode output looks, i.e. a series 
of A and B values alternating at half the frequency of the phase 
switch; in the big graph, instead, the measured values relative 
to the two switch states are shown separately for better clarity. 

From the ratio of the average levels of the two data streams 
we have calculated a value r - 0.96075 that has been applied 
to produce the differenced noise stream shown in the top panel 
of the right part of Fig.0 The graph shows that the offset has 
been removed as well as the long-timescale instabilities that 
are evident in the total power data. The power spectrum of the 
differenced data shows a final knee frequency of the order of ~ 
40 mHzl 

We also calculated r by taking the ratio of the standard 
deviation of the data streams, and by directly minimising the 
final knee frequency. In the first case we obtained a value of 
r = 0.95568, while in the second case we found r - 0.96076 
which is practically coincident with the one calculated with the 
ratio of signal levels. 

In Fig|S]we show the behaviour of the knee frequency ver- 
sus r around the minimum; the continuous line represents a 
polynomial fit used to determine the minimum. 



the knee frequency was defined as the frequency at which the 
power spectrum has a value equal to V2 of the average white noise 
level 



signal (uncallbrated units) Signal (uncallbrated units) 



Difference, A - r * B 




200 400 600 800 1000 0,001 0.01 0.100 1.000 10.000 100.000 

time (s) frequency (Hz) 

Fig. 7. Measured data from the 30 GHz Planck-LFI prototype 
radiometer The graph on the left shows about 15 min of to- 
tal power noise data from one of the two radiometer channels. 
The small inset shows how measurements actually appear at the 
output of each detector, i.e. as an alternating series of A and B 
values. The two panels on the right show the differenced noise 
stream in time and frequency domain. The final knee frequency 
is -40 mHz. 

4 (Hz) 

0.0415 1 ' ' ' ' \ ' ' ' ' \ ' ' ' ' \ ' ' ' ' \ ' ' ' ' \ ' ' ' ' 




0.0395 I : : : : \ : : : : \ : : : : \ : : : : \ : : : : \ : : : : 

0.9595 0.9600 0.9505 0.9610 0.9615 0.9520 0.9525 

r 

Fig. 8. Behaviour of /k versus r around the minimum. The 
curve was obtained by applying different values of r to the total 
power data shown in Fig. and then calculating the knee fre- 
quency of the differenced data. The continuous line is a poly- 
nomial fit of the discrete data points. 

Although the knee frequency obtained from the differential 
data is within the Planck-LFI requirements its value is higher 
than the theoretical value (< 10 mHz) that can be estimated 
from Eqs. and 

Part of this discrepancy was caused by thermal drifts in the 
300 K environment (that were observed to be of the order of 1 .3 
K/hr and can be clearly seen in the total power data streams); to 
evaluate this effect we removed the linear trend from the data 
and then calculated the spectrum of the resulting differenced 
time stream. In Fig.|5]we show the total power and the differ- 
ential radiometric outputs after the removal of the linear trends; 
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the noise power spectrum shows that the final knee frequency 
was reduced from ~ 40 to ~ 26 mHz. 

Signal (uncallbrated units) Signal (uncallbrated units) 



Difference, A - r * B 




200 400 600 800 1000 0.001 0.010 0.100 1.000 10.000 100.000 



time (s) frequency (Hz) 

Fig. 9. Left panel: total power radiometer data after the re- 
moval of the linear trend. Right panel: differetial output in time 
(top) and frequency domains (power spectrum, bottom). The 
knee frequency was estimated ~ 26 mHz. 

The remaining discrepancy of about ~15 mHz was likely 
to be caused by other non-idealities in the experimental setup 
which have not been completely understood yet. 

These results demonstrate how the pseudo-correlation 
scheme adopted for the LFl radiometers is effective in reduc- 
ing the effect of amplifier instability to very low levels also in 
experimental conditions far from the nominal flight operations. 
Furthermore, we have shown that the gain modulation factor 
can be determined with very simple analysis of a limited por- 
tion of the total power radiometer data. 

6. Impact of main systematic effects on the 
calculation of r 

In this section we present a brief analysis of the impact of some 
systematic effects on the calculation of the gain modulation fac- 
tor. We can identify the following three possible classes of ef- 
fects that can impact on the calculation of r: 

1 . effects that induce a variation of the total power signal lev- 
els, which have an influence especially in the case r is cal- 
culated using the ratio of the average total power sky and 
load levels; 

2. effects that induce a variation of the standard deviation of 
the total power noise streams, which have an effect in the 
case r is calculated using the ratio of standard deviations; 

3. effects that change the low frequency spectral behaviour 
of the noise streams and limit the ability to derive r from 
the analysis of the differenced noise stream frequency spec- 
trum. 

Considering that the most accurate method to derive r is 
from the ratio of the average signal levels, we focus on the first 



class of effects, which can be subdivided into the two following 
sub-classes: 

1. effects that induce roughly the same variation in both the 
sky and reference load total power signals (symmetric ef- 
fects), and 

2. effects that introduce variations in only one of the two sig- 
nals (sky or reference load, asymmetric effects). 

Effects that belong to the first category are, for example, 
fluctuations induced by radiometer temperature instabilities 
and by bias voltage variations. Effects belonging to the second 
category are, for example, the CMB dipole, telescope tempera- 
ture fluctuations, reference load variations. 

6.1. Symmetric effects 

Let us consider a systematic effect that causes a variation 6T 
in both the sky and reference load total power signals. Clearly 
this is completely equivalent to the case of a constant offset 
introduced by the back-end electronics discussed in Sect. 14.11 
Therefore we can conclude that we need systematic effects with 
an amplitude of the order of 1 K (which is 2-3 order of magni- 
tudes higher compared to what is expected for Planck-LFI) or 
more to to introduce uncertainties in the value of r greater than 
1%. 

6.2. Asymmetric effects 

Let us now consider a systematic effect that causes a variation 
(6Ty^^ only in the sky power signal. The relative uncertainty 
on the value of r determined by this effect is: 

6r (dTY^y 

— - ■ (26) 

'"o ^sky + Ta 

In the case we have a systematic variation only in the refer- 
ence load signal then the uncertainty on r is: 

6r {STy^ 

— = ^ — 7 (27) 

In Table m we report, for the Planck-LFI 30 and 100 GHz 
channels, the values of the sky and reference load variations 
that would determine a change in r equal to the needed accu- 
racy (first two rows), together with estimates of the variation 
in the sky and reference load signals due to the major expected 
asymmetric effects (last three rows). 

The values reported in Table|3show that in the Planck-LFI 
case the calculation of the gain modulation factor from the ratio 
of the average signal level is largely immune from any expected 
systematic variations in the sky signal and/or in the reference 
load. This means that in this case r can be considered constant 
at first order and the only significant changes are expected from 
long timescale variations of the radiometer noise properties. 
These changes will be recognised in the ground data analysis 
and compensated for by periodically updating the value of r 
during data reduction. 
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30 GHz 


100 GHz 




values in mK 




±146 


± 112 




±154 


± 113 




±3.4 


± 2.7 


(sn,i 


< ±1 


< ±1 




< ±1 


< ±1 



Table 2. Impact of asymmetric effects on the calculation of 

the gain modulation factor for the Planck-LFI radiometers. The 
first two rows show the variation in the sky and reference signal 
that would cause a change in r equal to the needed accuracy of 
+ 1% at 30 GHz and of ±0.2% at 100 GHz. The last three rows 
list the expected sky and reference signal variations caused by 
the major sources of asynmietric effects (CMB dipole, tele- 
scope temperature fluctuations, reference load instabiUties). 



7. Conclusions 

High sensitivity CMB radiometric measurements require a 
very low susceptibility of the receiver to amphfier 1// noise, 
which can be obtained only by differential measurements. The 
pseudo-correlation differential radiometer is a receiver partic- 
ularly suitable for CMB radiometric measurements, allowing 
differential measurements without the need of an active front- 
end switch. 

In this scheme the best 1// noise suppression is reached 
when the radiometer is completely balanced. For slightly un- 
balanced configurations, e.g. when a cryogenic internal load is 
used as a reference signal, the effects introduced by the offset 
between the two inputs can be minimised after detection either 
by a variable DC gain stage that maintains the differenced out- 
put as close as possible to zero, or by multiplying in software 
one of the two signals by a factor, r, that must be properly tuned 
to balance the receiver. At first order the parameter r is equal 
to the ratio between the sky and reference "total power" signal 
levels, i.e. r x (fs^y + Tn)/(Tisf + Tn). With this value it is pos- 
sible to obtain differenced noise streams with a knee frequency 
of the order of few tens of mHz and make the radiometer in- 
sensitive to fluctuations in the back-end amplifier gain. 

The accuracy required in the calculation of r from ground 
data must be determined by the instrument requirements on 1// 
noise and other systematic effects. In the case of Planck-LFI 
the main driver is the susceptibihty to thermal fluctuations in 
the warm back-end stage, which can be maintained to required 
levels if r is calculated with an accuracy ranging from +1% for 
the 30 GHz channel to ±0.3% for the 70 GHz channel. 

In our study we have considered three methods for calcu- 
lating r, using total power radiometer data. 

The most straightforward scheme uses the ratio of the aver- 
age level of the total power signals. This method proves quite 
simple, accurate and relatively immune from systematic effects 
like l/f amphfier fluctuations, thermal effects etc. In order to 
be applied it is necessary that any spurious offset introduced by 
the back-end electronics is known with a relative accuracy of 
the order of ±4%. This is the baseline method in the case of 
Planck-LFI. 



The second method uses the ratio of the standard deviations 
of the total power noise streams and is hmited by its sensitivity 
to l/f fluctuations of the front-end amplifier noise temperature, 
which sets a lower limit on the expected accuracy that is about 
one order of magnitude worse than method 1. 

The third method uses a minimisation strategy of the knee 
frequency of the differenced noise stream, by which it is possi- 
ble to obtain the same level of accuracy of method 1, provided 
that the lenght of the total power data stream is enough to re- 
solve the knee frequency of the differenced noise. 

The apphcation of this concept to prototype Planck-LFI 
radiometer data has shown that the "in software" offset balanc- 
ing method is effective and that the balancing condition can be 
calculated from a hmited amount of raw radiometric data in to- 
tal power mode. In particular by methods 1 and 3 applied to 
a 15 min data stream of experimental radiometric data it was 
possible to calculate the same (at the level of 0.001%) value of 
r. 

Although in Planck-LFI method 1 is considered as base- 
line, the other methods will be used during the ground data 
analysis for cross-checks, as they are sensitive to non ideahties 
in different ways. 

Further studies will be aimed at a better understanding of 
the impact of systematic effects on the gain modulation factor 
accuracy and to continue the analysis of laboratory radiome- 
ter data in order to check our predictions. In addition we will 
continue with more realistic and detailed simulations. 
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